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ABSTRACT 


In  this  report  the  fundamental  frequency  of  a  sandwich  composite  plate  containing 
woven  fabric  layers  is  determined  using  the  Rayleigh-Ritz  technique.  The  plate  contains  a 
Baltek  balsa  wood  core  and  face  sheets  made  of  both  unidirectional  and  woven  fiberite 
layers.  The  flexural  stiffnesses  for  the  face  sheets  are  obtained  using  the  mosaic  model.  Of 
specific  interest  is  the  fundamental  transverse  bending  mode  and  is  numerically  computed 
using  MATHEMATICA  for  simply  supported  boundary  conditions  and  clamped 
boundary  conditions  which  are  the  in  service  boundary  conditions.  Experimental 
frequencies  are  determined  for  simulated  free  boundary  conditions  and  are  used  as  a  means 
of  comparison  for  the  numerically  computed  fundamental  frequencies. 
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IV 


Introduction 


Sandwich  structures,  coniposite  sandwich  structures  in  particular,  have  been  used 
widely  in  the  many  industries  including  the  automotive,  space  and  aerospace  industries. 

As  a  result,  research  efforts  have  provided  an  understanding  of  the  both  the  buckling  and 
vibration  behavior  of  fundamental  structures  including  plates  and  shells.  Of  particular 
importance  is  the  bending,  buckling  and  vibration  of  composite  sandwich  plates.  Mirza 
and  Li  (1995),  for  example,  study  the  vibration  of  isotropic  sandwich  plates  by  employing 
Green’s  function  analysis.  Al-Qarra  (1989)  analyzes  the  deflection  of  both  sandwich 
beam  and  plates.  Finally,  Shi  and  Tong  (1994)  studies  the  localized  buckling  of  composite 
honeycomb  structures.  Now,  these  type  of  structures  are  becoming  popular  for  use  on 
both  naval  and  marine  structures  as  deck  plating,  bulkheads  and  columns.  In  this  paper, 
the  necessary  details  are  presented  for  the  Rayleigh-Ritz  analysis  which  can  be  used  as  a 
technique  to  numerically  determine  the  fundamental  frequency  of  composite  sandwich 
plates  containing  woven  laminas.  Experimental  results  are  presented  for  the  specific  plate 
described  herein. 
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Problem  Statement 


The  equations  governing  the  vibration  of  a  sandwich  plate  are  given  by  three 
partial  differential  equations.  For  the  composite  sandwich  plate  containing  orthotropic 
face  sheets  these  equations  are 
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where  Djj  represent  the  bending  stiffiiesses  with  D,j  =  Djj  +  Dgg,  G13  and  G23  are  the  shear 
moduli,  I  represents  the  rota^  inertia,  p  provides  the  density,  and  h  is  the  thickness  of  the 
core.  Also  appearing  in  Eq.  (1)  are  the  kinematic  variables  w,  ij;^  ,and  1);^  which 
correspond  to  a  first  order  shear  deformation  theory.  Appropriate  boundary  conditions 
must  be  appended  to  Eq.  (1)  in  order  to  have  a  well  posed  problem  and  can  be  either 
simply  supported,  clamped,  free  or  a  combination  thereof.  For  the  plate  consider  herein 
simply  supported  and  clamped  boundary  conditions  are  considered  and  are  given  by 

w  =  0  =  0  =  0  (2) 


for  simply  supported  and 

w  =  0  T'=0 

s  n 


(3) 


for  clamped  boundary  conditions. 
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Method  of  Analysis 


The  are  many  techniques  which  can  be  utilized  to  solve  the  equations  given  in  Eqn. 
(1).  Popular  approaches  include  finite  difference  and  finite  element  methods.  An 
approximate  approach  yielding  the  fundamental  frequencies  is  the  Rayleigh-Ritz  technique 
and  is  the  chosen  approach  for  this  problem.  To  this  end,  the  total  potential  energy 
corresponding  to  Eq.  (1)  is 

^  ^22^y,y  '^^\2!^x^y,y  ^(>e^^x,y  'Ky.x 

+  +  /j(Gi3W,/+G23W,/T  - 

The  displacement  parameters,  i|;,,  t|;y,  and  w  are  assumed  separable  and  can  be  expanded 
in  an  admissible  basis  given  as 

V  (=1  y=l 

%(x,y)  =  'tf:^fM)Dj(y)e‘'^‘  (5) 

1=1  M 

y<x,y)  =  EEY,y^iW^;CyV'“' 

i=l  ;=1 

Substituting  Eqn.  (5)  into  (4)  and  minimizing  provides  the  discrete  eigenvalue  problem  in 
the  form  of 

[K]{u}  =  0)2[M]{m}  (6) 

where  {u}  =  {c£y,py,Yij}^  is  an  ordered  triple  of  expansion  coefficients  corresponding  to  the 
eigenvector,  to  is  the  circular  natural  frequency,  and  [K]  and  [M]  are  the  stiffriess  and 
mass  matrices.  The  stiffness  matrix  [K]  is  of  the  form 
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and  mass  matrix  [M]  is  of  the  form 

[mV]  [m‘^]  [m^^l 

[M]  =  [m‘2]  [m22]  [^23]  (9) 

[m*^]  [m^^]  [m^^] 

and  the  elements  of  the  matrices  are  given  by 
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^rksl  ~ 

\ksl  = 

^r^l  =  9(E,,E,W^,F) 
m^^l  =  [0]  forp  *  q 


(10) 


Above,  the  ( ‘ )  indicates  derivative  and  ( •,  • )  denotes  a  Lj  integration.  It  is  important  to 
note  that  as  N,  the  number  terms  taken  in  the  displacement  expansion,  approaches  infinity 
the  solution  from  Eq.  (6)  approaches  that  of  Eq.  (1).  Moreover,  since  the  order  of  this 
system  is  3N^  x  3N^  ,  the  size  of  the  system  becomes  large  for  relatively  small  N.  Here,  N 
is  limited  to  seven  terms  in  both  expansions  resulting  in  a  maximum  size  of  147  x  147  for 
this  study. 
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Plate  Geometry 


The  plate  under  investigated  is  a  half  scale  model  with  dimension  of  20  ft  long  by  6 
ft  wide.  The  stacking  sequence  of  the  plate  is  given  by 
[woven/90°/90°/woven/90°/90°/core]s  and  is  shown  in  Figure  1. 


0.025in 
0.01  Sin 
0.013  in 
0.025  in 
0.013  in 

0.013  in 

1.3  in 

Centerline 


Figure  1.  Sandwich  plate  lamination  sequence. 


As  mentioned  earlier  and  after  specializing  the  stacking  sequence,  the  bending  stiffnesses 
appearing  in  Eq.  (1)  contain  contributions  from  both  unidirectional  layers  and  woven 
roving  fabric  layers.  For  an  arbitrary  sandwich  plate  with  identical  upper  and  lower  face 
sheet  construction  containing  N  layers,  Whitney  (1987)  indicates  the  face  sheets  are  the 
only  portion  of  the  plate  resisting  the  bending  stresses  and  hence  the  bending  stif&iesses 
can  be  computed  through 
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(11) 
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In  Eqn.  (1 1),  provides  the  coordinate  to  the  k  lamina  and  (2, /are  the  transformed 
stiffnesses  for  this  lamina  [4].  The  properties  for  each  material  are  given  in  Table  1. 


There  are  several  methods  which  can  be  used  to  model  the  flexural  stif&iesses  for 
the  woven  fabric  laminas.  Among  these  approaches  are  the  mosaic  approach  and  the  fiber 
undulation  approach.  The  mosaic  approach  models  the  woven  fabric  as  an  assemblage  of 
strips  in  orthogonal  directions  identified  as  the  warp  and  fill  directions  and  does  not  take 
into  account  the  undulating  nature  of  the  fibers  in  there  respective  directions  (Chou  1989). 
This  technique  effectively  models  the  woven  layer  as  a  cross-ply  laminate.  The  fiber 
undulation  model  does  take  into  account  the  undulating  nature  of  the  fibers.  The 
approach  produces  flexural  stiffness  constants  that  are  position  dependent.  Since  each 
face  sheet  has  a  thickness  much  smaller  than  thickness  of  the  core,  the  variation  in  height 
of  the  fibers  will  not  affect  the  overall  stiffness  of  the  face  sheets.  Therefore,  the  mosaic 
model  provides  a  reasonable  approximation  of  the  stiffnesses  for  sandwich  type  structures. 


Experimental  Configuration 

With  a  panel  of  this  size  and  weight,  it  was  impractical  to  construct  a  support  that 
could  approximate  a  good  fixed  boundary  on  aU  sides  of  the  panel.  With  the  identification 
of  the  fundamental  transverse  frequency  being  of  importance,  the  panel  was  suspended 
on  rubber  bungee  cords,  simulating  as  near  as  possible  free  boundaries.  The  bungee  cords 
were  passed  under  the  panel,  and  around  a  steel  test  frame.  A  total  of  22  vertical  cords 
was  used,  giving  a  load  in  each  bungee  in  the  range  20-30  (lbs).  Experience  dictates  that 
this  load  range  in  these  particular  bungee  cords  is  optimum  for  providing  good  isolation 


-7- 


and  very  low  suspension  natural  frequencies.  In  this  case  the  rigid  body  modes  were  aU 
below  about  1-2  Hz.  It  is  believed  that  these  modes  did  not  alter  the  analysis  to  any 
significant  extent.  A  schematic  of  the  tests  configuration  is  shown  in  Figure  2. 
Complete  details  of  the  experimental  investigation  can  be  found  in  Ratcliffe  (1996). 


Figure  2.  Schematic  -  Experimental 
Configuration. 

The  excitation  was  impact,  referenced  to  a  fixed  accelerometer.  Calibrated  FRF 
data  in  the  frequency  range  0-200  Hz  were  captured  for  impact  at  each  of  the  147 
coordinate  locations  on  the  panel.  A  main  consideration  during  testing  was  the 
accelerometer  location,  which  was  located  at  a  coordinate  that  was  close  to  one  edge  and 
well  away  from  the  center-line  and  lines  of  symmetry.  This  coordinate  was  unlikely  to  be 
on  a  node  line  at  resonance,  a  requirement  for  a  good  location.  While  there  may  be 
additional  natural  frequencies  that  were  not  identified,  the  results  of  the  analysis  (not  fully 
reported  here)  identified  aU  the  "uniform  plate"  modes  in  the  frequency  range  of  interest. 
Therefore,  it  is  believed  the  transducer  location  was  adequate  for  this  project. 
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Results  and  Discussion 


Nutrorical  results  from  the  Rayleigh-Ritz  analysis  have  been  generated  for  simply 
supported  (SSSS)  and  clamped  (CCCC)  boundary  conditions  on  all  sides.  The  Ritz 
technique  was  programmed  in  MATHEMATICA.  Convergence  results  are  only  necessary 
for  the  CCCC  since  the  exact  basis  functions  can  be  selected  for  the  SSSS  boundary 
condition.  Table  1.  shows  the  convergence  of  the  fundamental  transverse  bending  mode 
clamped  boundary  conditions.  The  fundamental  frequency  for  this  condition  at  N  =  7 
terms  is  7.839  Hz.  The  fundamental  frequency  for  the  simply  supported  boundary 
condition  is  4.152  Hz 


N 

(o(Hz) 

1 

13.146 

3 

8.035 

5 

7.881 

7 

7.839 

Table  1.  Convergence  results  for  the  transverse  frequency  for  clamped  boundary 

conditions  on  all  sides 

The  FRF  data  were  subject  to  a  modal  analysis  using  the  commercial  STAR 
program.  The  imaginary  component  of  the  FRF  is  normally  a  good  indicator  of  resonance, 
and  therefore  the  averaged  (imaginary)^  accelerance  was  initially  used  as  an  indicator  for 
natural  frequencies.  Several  frequency  bands  were  selected  for  detailed  analysis.  The 
panel  had  some  small,  resonant  fixtures  which  increased  the  number  of  resonances  above 
that  expected  for  a  uniform  panel.  In  all,  there  were  39  natural  frequencies  in  the  range  0- 
200  Hz.  The  lowest  ones  are  shown  in  Table  2,  which  includes  the  identified  modal 
viscous  damping  ratios. 
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Natural 

Frequency 

(Hz) 

Viscous 

Damping 

Ratio 

(%) 

7.55 

3.01 

11.87 

2.14 

20.42 

1.42 

25.76 

1.59 

Table  2.  Experimental  fundamental  frequencies  results  for  the  suspended  sandwich  plate 

The  mode  shapes  for  these  resonances  are  shown  in  Figures  3  and  4.  Fig  3.  shows  the 
fundamental  mode  shape  for  this  panel.  This  shape  is  identified  as  a  (2,  0)  mode,  since 
there  are  two  phase  changes  along  the  length  of  the  panel,  and  no  phase  change  across  the 
width.  Figs  4  show  the  next  mode  shape  and  is  identified  as  the  (1,1)  mode. 


Figure  3.  Graphical  presentation  of  fundamental  transverse  bending 

mode  of  7.6  Hz 
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Figure  4.  Graphical  presentation  of  fundamental  transverse  bending  mode  of 

11.9  Hz 


Comparison  of  the  numerically  computed  fundamental  frequency  with  the  experimentally 
determined  frequency  indicate  a  3.8  %  percent  difference  with  the  Ritz  frequency  being 
higher. 

Conclusions 

In  this  paper  the  fundamental  transverse  frequency  of  sandwich  plates  containing 
woven  fabric  layers  has  been  computed  using  the  Rayleigh-Ritz  method  and  through 
experimental  methods.  The  frequency  was  determined  for  both  simply  supported  and 
clamped  boundary  conditions  analytically  and  determined  experimentally  using  pseudo- 
free  boundary  conditions.  Modeling  the  woven  fabrics  using  the  mosaic  approach  is  a 
viable  method  for  sandwich  plates.  Comparison  between  the  two  methods  indicates  a 
good  agreement  with  only  3.8  %  difference  for  the  fundamental  transverse  bending  mode. 


Recommendations 


As  a  result  of  the  analysis,  the  following  recommendations  are  provided  for  future 
considerations: 

(a) .  Include  both  elastic  and  rotational  boundary  effects; 

(b) .  Optimization  for  minimum  fundamental  frequency  for  specified  geometry. 

(c) .  Investigation  full  scale  vibration  characteristics. 
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